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Amorphous solids are mechanically rigid while possessing a disordered structure similar to that of 
dense liquids. Recent research indicates that dynamical heterogeneity, spatio-temporal fluctuations 
in local dynamical behavior, might help understanding the statistical mechanics of glassy states. 



THE PUZZLE POSED BY AMORPHOUS 
MATERIALS 



From the point of view of statistical physics, glasses 
are mysterious materials. Glassy materials possess a me- 
chanical rigidity which is similar to the one of a crys- 
talline material. In a crystal, rigidity is a direct conse- 
quence of the long-range periodic order: it is not possi- 
ble to move a single particle in a perfect crystal (while 
preserving the crystalline order) without also moving an 
extensive set of neighbors, see Fig. HJa). While mechani- 
cally rigid, glasses do not seem characterized by any type 
of long-range order, see Fig. [ljb), they actually resemble 
ordinary dense liquids. The comparison between crystals 
and glasses suggests that perhaps a more subtle symme- 
try breaking takes place during the formation of a glass, 
one that is not obvious to the naked eye. This conun- 
drum has been a long-standing issue in condensed matter 
physics 

Experimentally, one faces the fundamental difficulty 
that liquids approaching the glass transition (for exam- 
ple by decreasing temperature) become too viscous to 
flow on experimental timescales, and fall out of thermal 
equilibrium without encountering any reproducible ther- 
modynamic phase transition. It is of course tempting to 
interpret this dramatic dynamic slowing clown as orig- 
inating from an underlying phase transition or critical 
point. Near an ordinary critical point, large scale spa- 
tial fluctuations develop, as the density fluctuations in 
the example shown in Fig. QJc), and the dynamics slows 
down Q. However, no such fluctuations are detected 
near the glass transition (see Fig.QJb)). Therefore, find- 
ing convincing evidence of an underlying phase transition 
governing the physics of amorphous materials would rep- 
resent an important progress in this field. 

In the last decade these questions have also come 
up in the field of soft condensed matter, in which dis- 
ordered structures known as "jammed" materials Q 
(foams, emulsions, colloidal gels, sandpiles) stop flow- 
ing when their density becomes large, without possessing 
long-range crystalline order, just like molecular glasses. 
J.D. Bernal M in the 60' s was one of the first physicists to 
suggest that disordered atomic fluids and granular pack- 
ings could be investigated using similar tools and, per- 
haps, understood using similar theoretical concepts-an 
idea that has since remained highly popular [f| . 



Two decades of research on dynamic heterogeneity (to 
be defined shortly) in amorphous materials have estab- 
lished that the formation of rigid amorphous structures is 
indeed accompanied by nontrivial spatio-temporal fluc- 
tuations, which become stronger as the glassy phase is 
approached and are characterized by growing dynamic 
correlation length scales Q. In this article we revisit 
the mounting evidence - using mostly the example of su- 
percooled liquids, where dynamic heterogeneity has been 
most widely analyzed - that the formation of amorphous 
materials is a complex collective phenomenon, which 
shares more similarities with ordinary critical points than 
the featureless structure shown in Fig. [TJb) would sug- 
gest. 



WHAT IS DYNAMIC HETEROGENEITY? 

The concept of dynamic heterogeneity as a key fea- 
ture that characterizes disordered materials has slowly 
emerged from experimental studies of highly viscous 
molecular liquids approaching the glass transition. In 
these systems, relaxation spectra measured through me- 
chanical or dielectric probes span a very broad range of 
relaxation times and are strongly nonexponential. This 
suggests the existence of wide distributions of relaxation 
rates. 

What is the microscopic origin of these broad distri- 
butions? Looking again at the disordered structure in 
Fig. \V[b), it is natural to imagine that the presence of 
structural disorder implies that atoms in different envi- 
ronments move differently. The physical picture is that, 
at any given time, different regions in a liquid might relax 
in a different manner and at a different rate, thus pro- 
ducing broad distributions of relaxation times. However, 
when the system is close to the glass transition but not 
yet a glass, particles constantly move and rearrange, and 
so the distinction between different spatial environments 
can only hold over a finite duration. 

Dynamic heterogeneity refers to the existence of tran- 
sient spatial fluctuations in the local dynamical behavior. 
Dynamic heterogeneity is observed in virtually all disor- 
dered systems with glassy dynamics Q ■ 

Direct confirmation of the heterogenous character of 
the dynamics stems for instance from computer simula- 
tions of simple models of supercooled liquids, see Fig. [2] 
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FIG. 1: (a) A periodic crystalline structure does not flow because preserving the crystalline order requires moving an extensive 
set of particles, (b) A mechanically rigid glassy structure exhibits neither the long-range order of a crystal nor the large scale 
density fluctuations observed at an ordinary critical point, (c) Large scale critical density fluctuations near the critical point. 



Whereas the averaged displacement of a particle in a 
given time window of width t is a smooth increasing func- 
tion of t, the time signals for individual particles shown 
in Fig. [5] have two important characteristics. 

(i) They are highly intermittent in time, being com- 
posed of a succession of long periods of time where par- 
ticles vibrate around well-defined locations, separated by 
rapid "jumps" . The waiting times separating successive 
jumps are statistically broadly distributed. 

(ii) The trajectories differ widely from one particle to 
another in the same system at the same time. Some 
particles undergo many jumps and move large distances 
while some other particles are nearly immobile over the 
entire time window. 

As all other features related to dynamic heterogene- 
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ity, such observations cannot be made from ensemble- 
averaged measurements. Indeed, the characterization of 
dynamic heterogeneity requires the development of ex- 
perimental techniques that are not only sensitive to aver- 
aged or typical behaviour, but can also resolve the fluc- 
tuations [3, Q. For theoreticians, the existence of dy- 
namic heterogeneity implies that fluctuations need to be 
taken into account in the description of transport proper- 
ties. Therefore, materials close to a glass transition differ 
qualitatively from ordinary fluids, where fluctuations can 
typically be neglected. 

The intermittency of single-particle trajectories, while 
a clear indication of spatio-temporal fluctuations, do not 
tell us how these fluctuations arc correlated in space. 
This point was first addressed in pioneering works us- 
ing four-dimensional NMR [i| and direct probing of fluc- 
tuations at the nanoscopic scale using atomic force mi- 
croscopy techniques in polymeric glasses [lfj| . Direct vi- 
sualizations of molecular trajectories are not yet possi- 
ble but recent, very elegant, experimental approaches 
using single molecule spectroscopy are getting close to 
it [ll|, |l2j. Direct visualisation is possible for differ- 
ent p types^ of glasses, such as colloidal [13| and granu- 
lar 



15| assemblies. All these spatially resolved mea- 
surement indicate that extended regions of space indeed 
transiently behave as fast and slow regions. 



NEED FOR HIGH-ORDER DYNAMIC 
CORRELATION FUNCTIONS 



FIG. 2: Time resolved squared displacements of individual 
particles in a simple model of a glass-forming liquid. The 
average behavior is shown as the smooth full line. Individual 
trajectories are composed of long periods of time during which 
particles vibrate around well-defined positions, separated by 
rapid jumps that are widely distributed in time and occur at 
different times and frequencies for different particles. 



We have described dynamic heterogeneity as spatial 
variations in the local relaxation rate. Although hard 
to detect in an experiment, these fluctuations are easily 
observed in a computer simulation, where the position 
of each particle is known at each time step. In Fig. [3] 
we show an example of the visualisation of the spatially 
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FIG. 3: Spatial map of single particle displacements in the 
simulation of a Lennard- Jones model of a supercooled liquid in 
two spartial dimensions. Arrows represent the displacement of 
each individual particle in a trajectory of duration comparable 
to the typical structural relaxation time. This map reveals 
that particles with different mobilities are spatially correlated. 

heterogeneous dynamics in a simple model of a super- 
cooled liquid, in two spatial dimensions. This visualiza- 
tion shows the existence of spatially extended "domains" 
where the amplitudes of single-particle displacements are 
correlated. Remember that these domains have no obvi- 
ous counterpart in the density fluctuations, Fig.QJb), and 
only appear when dynamics is considered. However, the 
spatial fluctuations in Fig. [3] arc obviously reminiscent of 
the critical fluctuations in Fig. [He), with one major dif- 
ference: while a thermodynamic quantity (the density) 
becomes critical in ordinary critical phenomena, fluctua- 
tions are only detected through dynamical quantities in 
highly viscous liquids. 

To characterize the spatial fluctuations shown in Fig. [3] 
in a statistical manner, one must resolve dynamics in 
both space and time and quantify deviations from the 
average behavior. To this end, one defines a "mo- 
bility", Ci(t, 0), which quantifies how much particle i 
moves between times t = and t [for instance Ci(t) = 
e -|ri(t)-rj(o)| Id ^ w jth d of the order of the particle size]. 
Given two particles at separation r, one can measure the 
degree to which their mobilities are correlated. It is con- 
venient to define a "mobility field" for a system composed 
of N particles via 

TV 

c(r;t > 0)=53cj(t,0)*[r-r i (0)] > (1) 

where Yi(t) is the position of particle i at time t. The 
spatial correlations of the mobility are finally captured 



by the correlation function 

G 4 (r; t) = (c(r; t, 0)c(0; t, 0)) - <c(r; t, 0)) 2 , (2) 

which depends only on the time t and the distance r = |r| 
as long as the ensemble average, denoted with brackets, 
is taken at equilibrium in a translationally invariant sys- 
tem; Gi(r;t) is known as a 'four-point' dynamic correla- 
tion function because it measures correlations of motion 
between time and time t arising at two points, and r. 

The analogy with fluctuations in critical systems be- 
comes clear in Eq. ([2]) if one considers the mobility field 
c(r; t, 0) as playing the role of the order parameter for the 
transition, characterised by nontrivial fluctuations and 
correlations near the glass transition. This analogy is 
now fully exploited in modern theoretical treatments (H . 

The above definition of a real-space correlation func- 
tion of the mobility represents a vital advance in the char- 
acterization of dynamical heterogeneity. For instance, 
it allows the language of field theory and critical phe- 
nomena to be used in studying dynamical fluctuations in 
glassy systems , [13] • Like in critical phenomena, if one 
makes the hypothesis that there exists a single dominant 
length scale £4 then one expects that for large distances 
the correlation function Gi{r\ t) decays as 

GA{r . t)K m e -r/Ut)^ (3) 

with p a critical exponent. It is also natural to define the 
susceptibility associated with the correlation function 

X4(t) = J d d vG 4 (r;t). (4) 

If the prefactor A(t) in Eq. ^ is known, the susceptibil- 
ity Xi(t) can be used to infer the typical number of par- 
ticles involved in the correlated motion shown in Fig. m 
That is, X4(t) ma y be interpreted as the "volume" of the 
correlated clusters. 

Further, Xi(f) can a l so be obtained from the fluc- 
tuations of the total mobility C(t, 0) = J d d r c(r; t, 0), 
through 

Xi (t) = N[(C(t,0f) - (C(t,0)f}. (5) 

In practice, this formula allows an efficient measure of the 
degree of dynamical heterogeneity, at least in computer 
simulations and in those experiments where the dynam- 
ics can be spatially and temporally resolved. As long 
as c(r; t, 0) appropriately quantifies atomic motion, x^{t) 
can be measured in a variety of systems, serving as a basis 
for fair comparisons of the extent of dynamical hetero- 
geneity and has become a central tool in characterizing 
dynamic heterogeneity in amorphous materials 0- 
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FIG. 4: Time dependence of the susceptibility X4-(fi) that 
quantifies spontaneous fluctuations of the relaxation function 
in a molecular dynamics simulation of a supercooled liquid. 
For each temperature, X*W nas a maximum, which shifts to 
larger times and has a larger value when T is decreased, re- 
vealing the increasing length scale of dynamic heterogeneity 
in supercooled liquids that approach the glass transition. 



FOUR-POINT SUSCEPTIBILITIES IN 
MOLECULAR, COLLOIDAL, AND GRANULAR 
GLASSES 



The dynamical function Xi{t) nas now been measured 
in computer simulations of many different glass-forming 
liquids, by molecular dynamics, Brownian, and Monte 
Carlo simulations [18|-|23j. An example is shown in Fig. [4] 
for a Lennard- Jones numerical model, but the qualitative 



behavior is similar in all cases 24-261: function of 
time Xi(~t) increases at first, it has a peak on a timescale 
of the order of the typical relaxation time of the fluid, and 
then it decreases at large times. This time dependence 
simply reflects the the transient nature of the dynamical 
heterogeneity. 

The peak value of Xi(t) approximately measures the 
volume over which structural relaxation processes are 
correlated. Therefore, the most important result ob- 
tained from data such as those presented in Fig. @] is 
the temperature evolution of the peak height, which is 
found to increase when the temperature decreases and 
the global dynamics slows down. Such data provides di- 
rect evidence that the approach to the glass transition 
is accompanied by the development of increasingly long- 
ranged spatial correlations of the dynamics. 

In experiments, direct measurements of X4(t) have 
been made in colloidal j2?| and granular materials |14l.[l5l| 
close to the colloidal and gra nular glass transitions, and 
also in foams [28| and gels |29j , because dynamics is more 
easily spatially and temporally resolved in those cases. 
The results obtained in all these cases are again broadly 
similar to those shown in Fig. HI both for the time de- 
pendence of Xii.^) an d its evolution with a change of the 
relevant variable controlling the dynamics. 



Obtaining information on the behavior of Xi{t) an d 
Gi(r;t) from experiments on molecular systems is diffi- 
cult, because it is hard to disentangle the spontaneous 
fluctuations embodied in Xi(t) m Eq. ([5]) from the ex- 
perimental noise. Such measurements are however cru- 
cial because numerical simulations and experiments on 
colloidal and granular systems can typically only be per- 
formed for relaxation times spanning at most 5-6 decades. 
On the other hand, in molecular liquids, up to 14 decades 
are in principle relevant, and extrapolation of simulation 
data all the way to the experimental glass transition is 
fraught with difficulty. Indirect estimates of Xi(t) from 
experiments are discussed below. 



REAL-SPACE MEASUREMENTS AND 
DYNAMIC STRUCTURE FACTORS 

We mentioned above that a growing peak in Xi{t) "di- 
rectly" reveals the growth of a dynamic correlation length 
scale as the glass transition is approached. This can only 
be true if the assumptions made in Eq. ([3]) for the scal- 
ing form of Gi{r\t) are correct. Dynamic length scales 
should in principle be obtained by direct measurements 
of a spatial correlation function. 

However, in contrast to X4(t), detailed measurements 
of G±{r;t) are technically more challenging as dynamic 
correlations must now be resolved in space over a large 
range of distances with a very high precision, and so there 
is much less data to draw on. From the point of view of 
numerical simulations where many measurements of X4 
were reported, the main limitation to properly measur- 
ing £4 is the system size. This might seem surprising as 
typical numbers extracted for the correlation length scale 
£4 are rather modest, but a precise determination of £4 
requires an accurate study of the tail of G±{r; t) at large 
r, which entails an important numerical effort [3oLl3ll|. 



However, such studies are important in that they allow 
the dynamical length scale £4^) to be measured directly. 
Moreover, such studies help infer the behavior of £4^) 
from measurements of Xii^)- Published work is consis- 
tent with X4(t)/G4(0, t) representing the number of par- 
ticles involved in heterogeneous relaxation. Therefore, 
truly 'direct' measurements indeed confirm that the in- 
crease of the peak of X4(*) corresponds, as expected, to 
a growing dynamic length scale £4^) [U [23. 26l 3fj|. 31 



Instead of direct inspection of G^r; t), it is often con- 
venient to analyse its Fourier transform, 



S 4 (q;t) = J d d re l 



r G 4 (r;i). 



(6) 



which is known as the four-point structure factor of dy- 
namic heterogeneity. In Fourier space, the large domains 
observed in Fig. [3] impact the low wavevector behavior of 
S^(q;t) in the form of a peak that grows when the glass 
transition is approached. This peak is often fitted with 
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the Ornstein-Zernikc functional form which is frequently 
used in conventional critical phenomena 0] . 



EXPERIMENTAL ESTIMATES OF MULTIPOINT 
SUSCEPTIBILITIES 

Although readily accessible in numerical simulations, 
the fluctuations that give access to Xi(t) are in general 
very small and impossible to measure directly in experi- 
ments, except when the range of the dynamic correlation 
is macroscopic, as in granular materials (l4| or in soft 
glassy materials where it can reach the micrometer and 
even millimetre range (2|| (2i|, 32 1. To access Xi(~t) m 



molecular liquids, one should perform time-resolved dy- 
namic measurements probing very small volumes, with a 
linear size of the order of a few nanometers. 

Fortunately, simpler alternative procedures exist. The 
central idea underpinning these solutions is the realiza- 
tions that if it is generally hard to detect noise in an 
experiment, it is usually simpler to measure the response 
of a system to an external perturbation. In the linear 
response regime, both types of measurements can often 
be related to one another by fluctuation-dissipation theo- 
rems Q . The physical motivation is that while four-point 
correlations offer a direct probe of the dynamic hetero- 
geneities, other multi-point correlation functions might 
also give useful information about the microscopic mech- 
anisms leading to these heterogeneities. 

For example, one might expect that a local fluctua- 
tion of the enthalpy 5h x (t = 0) at position x and time 
t = triggers or eases the dynamics in its surroundings, 
leading to a systematic correlation between Sh x (t = 0) 
and c(x +r;t, 0). A similar effect is expected for a lo- 
cal fluctuation of the density. These physical intuitions 
suggest the definition of a family of three-point correla- 
tion functions that relate thermodynamic fluctuations at 
one point to dynamical ones at another point. Crucially, 
some of these three-point correlations are both experi- 
mentally accessible and give bounds o r ap proximations 
to the four-point dynamic correlations 



Based on this insight, one may obtain a lower bound on 
Xi{t) by using a fluctuation-dissipation relation, which is 
valid at e quili brium when the energy is conserved by the 
dynamics [33j: 



cp 



IxAtf 



(7) 



where Xr(t) quantifies the response of the average mobil- 
ity to an infinitesimal change in the temperature T, and 
cp is the specific heat per particle. The response xr(i) 
can be experimentally accessed by measuring the average 
of a dynamical correlator, (C(f, 0))t, at nearby temper- 
atures, T and T + ST, in the linear regime ST <C T: 

(C(t,Q)) T+ST - (C(t,0)) T 



10 



10 1 



2 r 




x 10° 



io- 



-1 _ 



IO" 2 
10 



Silica 
Lennard-Jones 
Hard spheres - -X - - 
BPM ---A-- 
Glycerol - - -a - - 
OTP 
Salol 

Propylene carbonate - - <0- - 
3- Fluor o aniline 
Propylene glycol - 

B 2 3 - Wr- 

m-toluidine - - © - - 
Decaline ---€>-- 



io 6 

Ta/to 



10 



10 



10 



14 



XTit) 



ST 



(8) 



FIG. 5: "Dynamic scaling" relation between the number of 
dynamically correlated particles, evaluated by the peak height 
of X4, and relaxation timescale, r a , for a number of glass- 
formers [2f|, determined using the bound provided by Eq. (0. 
For all systems, dynamic correlations increase when the glass 
transition is approached. The full line through the data [3f| 
suggests a crossover from algebraic, \i ~ T ay to logarithmic, 
X4 ~ exp(T„), growth of dynamic correlations with increasing 



The main experimental advantage of ([8j is that spatio- 
temporal resolution is not needed, contrary to Eq. ([5]) . 

Detailed numerical simulations and theoretical argu- 
ments strongly suggest that the right hand side of 
Eq. actually provides a good estimate of %4 (t) in su- 
percooled liquids, and not just a lower bound. Similar 
estimates exist considering density instead of the tem- 
perature in Eq. ((SJ. These are useful when considering 
colloidal or granular materials where the glass transition 
is mostly controlled by the packing fraction. The quality 
of the corresponding lower bound was tested experimen- 
tally on granular packings close to the jamming transition 
34 1 , and numerically for colloidal hard spheres [10, |35[ . 

We show in Fig. [5] a compilation of data [3(§ for the 
evolution of the peak height of x&i in a representa- 
tion inspired by the theory of dynamic critical phenom- 
ena These data represent an experimental confir- 
mation that dynamic fluctuations and correlation length 
scales grow appreciably when molecular liquids approach 
the glass transition. However, we also learn that typ- 
ical length scales do not become very large (remaining 
in the nanometer scale for molecular glass-formers) be- 
fore liquids vitrify in a nonergodic state, and that an 
"ideal" glass critical point is not readily accessible to ex- 
periments. 
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THEORETICAL DEVELOPMENTS 

The above results are also relevant because many the- 
ories of the glass transition have assumed or predicted 
that the dynamics slows down because there are increas- 
ingly large regions over which particles have to relax in a 
correlated manner However, in the absence of exper- 
imental signs of growing length scales, these theoretical 
constructions would have remained speculative. 

The measurements of the spatial extent of dynamic 
heterogeneity, in particular Xi(t) an d G4(r;i), seem to 
provide the long-sought evidence that the glass transi- 
tion must indeed be considered as a form of critical or 
collective phenomenon involving growing timescales and 
lengthscalcs. This is important for the field of glass tran- 
sition, even though a conclusive understanding of the 
relationship between dynamical length scales and relax- 
ation timescales is still the focus of intense research Q ■ 

From a theoretical perspective, we are familiar with the 
idea, borrowed from equilibrium critical phenomena Q, 
that when correlation length scales get large, microscopic 
features of the system become unimportant and "univer- 
sal" behaviors emerge. Whether realistic glassy systems 
have length scales that are large enough for such a univer- 
sal description remains unclear. Although most theoret- 
ical approaches are in this spirit, one should perform an 
equally careful treatment of pre-asymptotic effects, which 
obviously matter for experiments that are performed far 
away from (putative) criticality. Therefore, theories of 
the glass transition are still crude descriptions of reality, 
despite large research efforts. 

Distinct microscopic mechanisms have been proposed 
that all give rise to growing dynamic correlations similar 
to the ones revealed by four-point functions. This rela- 
tive abundance of approaches can be interpreted in two 
ways. A first view is to consider dynamic heterogene- 
ity as a unifying physical phenomenon that accompanies 
the formation of glasses, and is central to its theoretical 
understanding. The second, more pessimistic interpreta- 
tion, is that if any reasonable theory predicts increas- 
ingly spatially heterogeneous dynamics, then studying 
dynamic heterogeneity might not help make substantial 
theoretical progress. This argument is balanced by the 
fact that the mere observation of dynamic heterogeneity 
has forced theoreticians to work hard on many fronts to 
understand what sort of microsopic mechanism is able 
(or not) to provide for instance new insights or quantita- 
tive predictions for the behavior of multipoint correlation 
functions, susceptibilities and dynamic correlation length 
scales 0- 

Several theoretical schemes, such as mode-coupling 
theory [37J , the random first-order transition theory (38| , 
dynamic facilitation approach [3911 . and frustration- 
limited domain scaling picture [40[ have now been de- 
veloped to the point that they provide useful theoretical 



guides to understand and analyze various aspects of dy- 
namic heterogeneity, see for a general theoretical re- 
view. One can thus hope that these approaches contain 
useful seeds for the construction of a "unified" theory 
of the glassy state. In all of these approches, dynamic 
heterogeneity and spatio-temporal fluctuations feature as 
central concepts. 



PERSPECTIVES 

Although we mostly discussed supercooled liquids near 
the molecular glass transition, we must emphasize that 
other materials have played an important role in develop- 
ing the concepts and tools described above. For instance, 
several soft condensed materials such as colloidal assem- 
blies have been instrumental in understanding the phe- 
nomenon of dynamic heterogeneity, because they share 
important similarities with molecular liquids but their 
elementary constituents arc so much larger than atoms 
(in the range 50 nm - l^m) that they can more easily be 
visualized |13ll27l |. Driven granular media (sheared or ag- 
itated systems) also undergo a "granular" glass transition 
which is empirically similar to the molecular transition. 
In that case, grains can even be directly tracked using a 
standard camera [lj,|l5|. These studies have contributed 
to give flesh to the concept of dynamic heterogeneity, be- 
cause spatial fluctuations could be seen. 

In the introduction we alluded to the rigidity transition 
occurring in athcrmal disordered granular packings. This 
"jamming" transition arises in an assembly of rigid par- 
ticles when the system cannot be compressed anymore, 
and is thus mainly a geometric transition where thermal 
fluctuations play no role. For spherical particles of equal 
sizes, this transition occurs near 'random close packing' 
<^ rcp « 0.64 The jamming transition is relevant also 
for athermal assemblies of soft particles, such as foams 
and emulsions, which are thus additional examples of dis- 
ordered rigid materials [3] • Connections with the physics 
of glasses are still rather speculative but are currently the 
focus on an important research effort. Detailed studies of 
dynamic heterogeneity in packings of soft and hard par- 
ticles near random close packing started to appear only 
recently [Il-44|, and could help elucidating similarities 
and differences between glasses and granular materials. 

One of the most frequently asked questions in studies 
of dynamical heterogeneity is whether the observed dy- 
namic fluctuations and correlations might have a struc- 
tural origin: Is there, after all, a "hidden" thermody- 
namic order parameter which would exhibit spatial fluc- 
tuations comparable to the ones revealed by dynamic 
heterogeneity studies? This question has attracted sus- 
tained interest. For example, in very early numerical 
work on dynamic heterogeneity, immobile regions were 
discussed in connection with compositional fluctuations 
in fluid mixtures [45| . However, dynamic heterogeneity 
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would not have emerged as an important concept if a 
simple, direct connection between structural order and 
relaxation dynamics had been satisfactorily established 
in amorphous materials. In that case, research would be 
dedicated to understanding the development of structural 
correlations at low temperatures in supercooled liquids, 
and to developing tools to measure, quantify and analyse 
such static features. 

Having said this, recent research on isoconfigurational 
ensembles fjif. amorphous order [13] and point-to-sct 
correlations [48[ , suggests that the structure of disordered 
materials might well be the next topic where new discov- 
eries and concepts will emerge in the near future. While 
two-point static correlations are poorly correlated to the 
evolution of the glassy dynamics, there is plenty of room 
for inventing more complicated correlation functions that 
could more accurately characterize the local structure of 
complex disordered media and explain their fascinating 
physical properties. 
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